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We investigate BBN in scalar-tensor theories of gravity with arbitrary matter couplings and self- 
interaction potentials. We first consider the case of a massless dilaton with a quadratic coupling to 
matter. We perform a full numerical integration of the evolution of the scalar field and compute 
the resulting light element abundances. We demonstrate in detail the importance of particle mass 
thresholds on the evolution of the scalar field in a radiation dominated universe. We also consider the 
simplest extension of this model including a cosmological constant in either the Jordan or Einstein 
frame. 
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I. INTRODUCTION 

The concordance model of cosmology calls for the 
introduction of a cosmological constant or a dark en- 
ergy sector. Various candidates have been proposed 0, 
among which the possibility that gravity is not described 
by general relativity on large cosmological scales. It is 
of interest therefore, to test our theory of gravity in a 
cosmological context. This can be achieved in two com- 
plementary ways, either by designing model independent 
tests (see e.g. Refs. 0EJ|j for a discussion of the various 
possible tests) or by considering a class of well motivated 
theories and use all available data to determine how close 
to general relativity we must be. 

Among all extensions of general relativity, scalar- 
tensor theories are probably the simplest in the sense 
that they consider only the introduction of one |5| (or 
many scalar field (s) universally coupled to mat- 
ter. These theories involve two free functions describing 
the coupling of the scalar field to matter and its self- 
interaction potential. They respect local Lorentz invari- 
ance and the universality of free fall of laboratory-size 
bodies. They are motivated by high-energy theories try- 
ing to unify gravity with other interactions which gener- 
ically involve a scalar field in the gravitational sector. 
In particular, in superstring theories the supermulti- 
plet of the 10-dimensional graviton contains a scalar field, 
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the dilaton, and other scalar fields, moduli, appear dur- 
ing Kaluza-Klein dimensional reduction of higher dimen- 
sional theories to our usual four dimensional spacctimc. 

In cosmology, two main properties make these theo- 
ries appealing. First, an attraction mechanism toward 
general relativity 0, Q was exhibited. This implies that 
even if the tests of general relativity in the Solar sys- 
tem set strong constraints on these theories, they may 
differ significantly from general relativity at high red- 
shift. Second, it was shown that the general mechanism 
of quintessence was conserved ^| |n| if the quintessence 
field was non-minimally coupled and that the attraction 
mechanism toward general relativity still held with run- 
away potentials p^ . ll3LIT3| . These extended quintessence 
models are the simplest theories in which there is a long 
range modification of gravity, since the quintessence field 
is light, and they allow for a very interesting phenomenol- 
ogy 

Cosmological data give access to various aspects of 
these models. The cosmic microwave background (CMB) 
tests the theory in the linear regime 0, 0, 0, 0, |2(| 
while weak lensing opens a complementary window on 
the non-linear regime [2lL |22| . Solar system experiments 
give information on the theory today and big-bang 
nucleosynthesis (BBN) allows us to constrain the at- 
traction mechanism toward general relativity || at very 
high redshift. 

BBN is one of the most sensitive available probes of 
the very early Universe and of physics beyond the stan- 
dard model. Its success rests on the concordance be- 
tween the observational determinations of the light ele- 
ment abundances of D, 3 He, 4 He, and 7 Li, and their the- 
oretically predicted abundances [U [24|. Furthermore, 
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measurements of the CMB anisotropies by WMAP |25j 
have led to precision determinations of the baryon density 
or equivalently the baryon-to-photon ratio, r\. As ij is the 
sole parameter of the standard model of BBN, it is possi- 
ble to make very accurate predictions |26ll27ll28l l29. 30] 
and hence further constrain physics beyond the standard 
model [all- 
in particular, the 4 He abundance is often used as a 
sensitive probe of new physics. This is due to the fact 
that nearly all available neutrons at the time of BBN 
end up in He and the neutron-to-proton ratio is very 
sensitive to the competition between the weak interaction 
rate and the expansion rate. Of interest to us here is the 
effect of modifications to gravity which will directly affect 
the expansion rate of the Universe through a modified 
Friedmann equation. 

The WMAP best fit assuming a varying spectral index 
is ttbh 2 = 0.0224 ±0.0009 and is equivalent to 7?io,cmb = 
6.14 ± 0.25, where 7710 = 10 10 t]. Using the WMAP data 
to fix the baryon density, the light element abundances 
HE HE HE HE HE can be quite accurately predicted. 
Some BBN results are displayed in Table 1. 

The effect of scalar-tensor theories of gravity on the 
production of light elements has been investigated ex- 
tensively (see e.g. Ref. HE for a review). As a first 
step, it is useful to consider only the speed up factor, 
£ = H/Hgr,, that arises from the modification of the 
value of the gravitational constant during BBN HE HE . 
Other approaches considered the full dynamics of the 
problem but restricted themselves to the particular class 
of Jordan-Fierz-Brans-Dicke theo ry 133 . of a massless 
dilaton with a quadratic coupling |35ll36j | or to a general 
massless dilaton HE- It should be noted that a com- 
bined analysis of BBN and CMB data was investigated in 
Ref. HE and Ref. HE • The former considered G constant 
during BBN while the latter focused on a non-minimally 
quadratic coupling and a runaway potential. We stress 
that the dynamics of the field can modify CMB results 
so that one needs to be careful while inferring fib from 
WMAP. 

The goal of this article is to implement scalar-tensor 
theories in an up-to-date BBN code. This will comple- 
ment our existing set of tools which allows us to con- 
front scalar-tensor theories with observations of type la 
supernovae and CMB anisotropies as well as weak 
lensing |2l| . In particular, the predictions to be com- 
pared with observations can be computed in the same 
framework for any self-interaction potential and matter- 
coupling function. 

We first recall, in § [nj the equations describing the 
theory to be implemented in our BBN code and we also 
discuss local constraints. As a check of our code, we 
consider, in § IIIII the case of a massless dilaton with 
quadratic coupling HE- ^ n particular, we perform a full 
numerical integration up to the present that can be com- 
pared with the analytical results of Ref. HE • We update 
the constraints on this model by taking into account the 
latest BBN data discussed above. We reaffirm that only 



helium-4 is sensitive to the modification of gravity con- 
sidered here. In § I1VI we will consider the simplest exten- 
sion of this model by introducing a cosmological constant. 
Such a constant can be introduced as a constant poten- 
tial either in the Einstein frame, hence keeping the dila- 
ton massless, or in the Jordan frame, hence generalizing 
the constant energy density component. Both cases are 
considered and we conclude in Section Applications 
to various cases of cosmological interest will be presented 
in a follow-up article. 



II. IMPLEMENTING SCALAR-TENSOR 
THEORIES OF GRAVITY IN A BBN CODE 

A. Scalar-tensor theories in brief 

In scalar-tensor theories of gravity, gravity is mediated 
not only by a spin-2 graviton but also by a spin-0 scalar 
field that couples universally to matter fields. In the 
Jordan frame, the action of the theory takes the form 



S 



d 4 x 



16ttG* 
+S m [g,j. v ; tp 



g[F{ V )R-g^Z{^^ th 



2t%)] 
(1) 



where G* is the bare gravitational constant from which 
we define = 8irG*. This action involves three arbi- 
trary functions (F, Z and U) but only two are physical 
since there is still the possibility to redefine the scalar 
field. F needs to be positive to ensure that the graviton 
carries positive energy. S m is the action of the matter 
fields that are coupled minimally to the metric g^ with 
signature (-, +, +, +). 

The action JQ| can be rewritten in the Einstein frame 
by performing the conformal transformation 



as 



S 



9l v = Hv)9nv 



-g* [R* - 2g% v dpiptdvip* 



d 4 x 
16ttG» 



+S m [A 2 (cp„)g* u ;ip] 



(2) 
(3) 



The field tp* and the two functions A(tp*) and V(tp*) are 
defined by 



dtp* 
dip 

A(<p*) 
2V(<p*) 



dlnFO) 
dp 



-1 2 



Zjp) 
2F{tp) 



(4) 

(5) 
(6) 



We will denote any Einstein frame quantities by a star 
(*), e.g. R* is the Ricci scalar of the metric g*^ v . The 
strength of the coupling of the scalar field to the matter 
fields is characterized by 



ot{ip*) 



din A 

dip* 



(7) 
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TABLE I: BBN results for the light element abundances assuming the WMAP-inferred baryon density. 



Source 


Y P 


D/H 


^He/H 


'Li/H 






xlO -5 


xl0~ 5 


xlO- 10 


Coc et al. (2004) 


0.2479±0.0004 


2.60 ±0.17 


1.04 ±0.04 


4.15 ±0.46 


Cyburt et al. (2003) 


U.44B4_ 0005 


z - 1 ^-0.16 




o 76 +i. tin 



and wc also define 



and conservation equations are given by 



da 

dip* 



(8) 



It is useful to study both the Einstein and Jordan 
frames. In the Jordan frame, matter is universally cou- 
pled to the metric. The Jordan metric defines the length 
and time as measured by laboratory apparatus so that all 
observations (time, redshift,...) have their standard inter- 
pretation in this frame. However, to discuss the theory it 
is often better to use the Einstein frame in which the ki- 
netic terms have been diagonalized so that the spin-2 and 
spin-0 degrees of freedom of the theory are perturbations 
of <?*„ and ip* respectively. The physical properties of 
both frames are of course identical. For example, when 
we refer to the time variation of the gravitational con- 
stant (in the Jordan frame), we have assumed fixed par- 
ticle masses. In contrast, in the Einstein frame, we would 
infer a fixed gravitational constant and varying masses. 
In both frames, the quantity Gm 2 (which is physically 
measureable) varies in the same way. 



Z(ip + 3Hip) = 3F^(H + 2H 2 ± 



K 
R 2 



--Z v p 2 - U v 



(12) 

p ± 3H (p ± P) = 0. (13) 
If we define the density parameters today by 

8irG*po 



n 



3H 2 F 



(14) 



the evolution of the energy density of a fluid with con- 
stant equation of state w = P/p takes the usual form 



87rG* 

where z is the redshift defined by 1 + z = Rq/R. 



2. Equations in Einstein frame 



(15) 



B. Friedmann equations 

1. Equations in Jordan frame 

We consider a Friedmann-Lemaitre universe with met- 
ric in the Jordan frame 



ds 2 = -di 2 ± R 2 (t)j ij dx i dx j 



(9) 



where 7^ is the spatial metric and R the scale factor. 
The matter fields are described by a collection of perfect 
fluids of energy density, p and pressure P. It follows that 
the Friedmann equations in Jordan frame take the form 



//•' • -^7 ) = 8irG*p+^Zp 2 -3HF + U(10) 
K 



-2F[H-^J = 8irG«(p + P) + Zip 2 

+F-HF 



(11) 



where a dot refers to a derivative with respect to the 
cosmic time t and H = dlni?/dt. The Klein-Gordon 



The scale factor and cosmic time in Einstein frame are 
related to the ones in Jordan frame by 

R = A(<p t )R t , dt = A(tp t )dt t (16) 
so that the redshifts are related by 



l + z= ^(l + z„). 
A 



(17) 



The Friedmann equations in this frame take the form 
3 (h 2 + J) = 87rG*p* ± V* + 2V{p,) (18) 



3 d 2 i?* 



4ttG*(p* ± 3P*) ± 2tJj 2 - 2V(ip*Xl9) 



Rl dtl 

where we have introduced H* = dlni?*/di* and 

%l)*=dip*/dU. (20) 

These equations take the same form as the standard 
Friedmann equations for a universe containing a perfect 
fluid and a scalar field. The Klein-Gordon equation takes 
the form 

^ ± 3#*V* — t~~~~ — 47rG,a(¥>.)(p, - 3P*) (21) 
dr* acp* 
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while the matter conservation equation is given by 
dp* 



dU 



+ 3H*(p* + -?*) = a{p*){p* - 3P*)V>*. (22) 



These equations differ from their standard form due to 
the coupling that appears in the r.h.s. The solution of the 
evolution equation i122[) can be obtained from the relation 
between the energy density and the pressure of a fluid in 
Einstein frame and their Jordan frame counterparts 



p* = A 4 p, P* = A 4 P 
which imply in particular, that 



(23) 



'•=i§K£r" + " ,(i+i - )s(,+ "' m 

for a fluid with a constant equation of state. 

C. Constraints today 

1. Post-newtonian constraints 

The post- Newtonian parameters (see Refs. can 
be expressed in terms of the values of a and (3 today as 



7 PPN -1 



2al 



1 + < 



2 • 







PPN 



1 = - 



1 (3 o>l 



2 (1 + <*§)*■ 



(25) 



Solar System experiments set strong limits on these pa- 
rameters. The perihelion shift of Mercury implies [4 



1 27 



PPN 







PPN 



1| < 3 x 1(T 3 , 



the Lunar Laser Ranging experiment \42\ sets 



47 



PPN oPPN 



-(0.7 ±1) x 10" 



(26) 



(27) 



Two experiments give a bound on 7 PPN alone, the Very 
Long Baseline Interferometer 



PPN 



II < 4 x 10" 



(28) 



and the measurement of the time delay variation to the 
Cassini spacecraft near Solar conjunction |44| 



7 PPN - 1 = (2.1 ±2.3) x 10- 5 . 



(29) 



These two last bounds imply ao to be very small, typi- 
cally a 2 , < 10~ 5 while Po can still be large 0]. Binary 
pulsar observations impose that (3q > —4.5. Note that 
even though (3q is not bounded above by experiment, we 
will assume that it is not very large, typically we as- 
sume Po ^ 100, so that the post-Newtonian approxima- 
tion scheme makes sense. 



2. Gravitational constant 

The Friedmann equations in the Jordan frame define 
an effective gravitational constant 



Gcff = G*/F — G*A 2 



(30) 



This constant, however, does not correspond to the grav- 
itational constant effectively measured in a Cavendish 
experiment. The constant measured in this type of ex- 
periment is 

G cav = G*A 2 {1 + a 2 ) (31) 

where the first term, G*Aq, corresponds to the exchange 
of a graviton while the second term, G^A^a 2 ,, is related 
to the long range scalar force. 

Assuming fixed particle masses, the time variation of 
the gravitational constant is bounded 01 by 



1 dG r 



|er | < 5.86 x lO -2 /! -1 . (32) 



Choosing the number of Einstein frame e-folds as a time 
variable, 



p = - ln(l + 2*), 



implies that 
2a 



00 



1 



00 
2 



dip* 
dp 



a . 



(33) 



(34) 



Note that the limit /3 = — (1 + a 2 ) that corresponds to the 
so-called Barker theory |48j in which A = cos ip* leads 
to a = whatever the value of a and ip'^ so that the 
gravitational constant is strictly constant even though 
gravity is not described by general relativity. 



D. Numerical implementation 

The nuclear reaction network takes its standard form 
in Jordan frame. To compute the light elements abun- 
dances during BBN, one only needs to know the expan- 
sion rate history, H(z), from deep in the radiation era 
up to today. It is thus convenient to express the Hubble 
parameter in the Jordan frame in terms of the one in the 
Einstein frame, using Eq. 1)16(1 . as 



AH = [H* + aip^tfj*} 



(35) 



where ip* is defined by Eq. 1(20(1 . Eq. 135(1 can also be 
expressed in the simple form 



AH = H* 



1 1 i \ dl P* 
1 + a{tp*) — 

dp 



(36) 
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It follows that, in terms of the cosmic time t, the equa- 
tions of evolution can be recast as 



d(p* 

~dT 

dtp* 

~~dT 



= A l (ip*)ip* 

= -A-^ip*) [ZH*ip* 

+ <iirG*a(ip*)A 4 (ip*) ^(1 - 3wi)pi 



dV 

d<p* 



(37) 



(38) 



hi = ^A^)j2Pi+yi+lvM-§m 



H = A^ 1 [H* + a{ip*)tp*) 



(40) 
(41) 



The numerical integration is performed as follows. 
First we choose some initial value ip- in *, Vw = deep 
in the radiation era (typically, z- m ~ 10 12 and we inte- 
grate the system (|37I41|I to z = 0. We perform a shoot- 
ing method so that the solution reaches the value Oao 
and Gat today, which fixes G* and the energy scale of 
the potential. At this stage the value of ip*o and «o are 
known. We also keep track of "00* to infer the time vari- 
ation of the gravitational constant, G cav , and check its 
compatibility with the constraint ll-il'l) . Subsequently, we 
perform a second integration of the same system includ- 
ing the nuclear reaction network. 



III. MASSLESS DILATON WITH QUADRATIC 
COUPLING 

The simplest model to consider consists of a massless 
dilaton with a quadratic coupling to matter. That is, 

1 



V(<p*)=0, 
It follows that 



.4 



a a(¥>.) 



= M- (42) 



ao = @<Po* 



(43) 



This model has been studied in detail in the literature, 
both in terms of its dynamics [1,0 and of its BBN predic- 
tions [33, Hg • We use it as a test model to check our nu- 
merical scheme. In particular, the analytical behaviour 
of the field during the radiation and matter eras after 
BBN was obtained for a flat universe without cosmo- 
logical constant in Ref. |35|. The numerical integration 
through BBN was also matched to this solution. Since 
we would like to use the same integration scheme for any 
potential and coupling, we can not rely on a particular 
analytic solution. It is used in this particular case only 
to check the accuracy of our code. 



A. General study 

As long as V = 0, the Klein-Gordon equation l|21|) can 
be rewritten in terms of the variable p defined by Eq. I|33|) 



3 — ip 



^'l + (1 - w)ti = -a(<p.)(l - 3w). (44) 



As emphasized in Ref. 8j, this is the equation of mo- 
tion of a point particle with a velocity dependent in- 
ertial mass, m(ip*) = 2/(3 — ip* 2 ), evolving in a po- 
tential a(ip*)(l — 3w) and subject to a damping force, 
— (1 — w)ip'*. During the cosmological evolution the field 
is driven toward the minimum of the coupling function. 
If > 0, it drives ip* toward 0, that is a — > 0, so that the 
scalar-tensor theory becomes closer and closer to general 
relativity. When (3 < 0, the theory is driven way from 
general relativity and is likely to be incompatible with 
local tests unless ip* was initially arbitrarily close to 0. 
Thus, we will restrict our analysis to f3 > 0. 

We need to consider three regimes : (i) deep in the ra- 
diation era, (ii) the effect of particle annihilation during 
the radiation era (electron-positron annihilation in par- 
ticular) and (iii) the transition between the radiation and 
matter era. 



1. Deep radiation era 

Deep in the radiation era, w = 1/3 and the coupling 
to ip* is not efficient. The equation of evolution reduces 
to 



3 - V? 

This can be integrated to give 



3^ 



0. 



<P* = ^p* 



V31n 



where cti is defined by 



<P*i 



3 - 



(45) 



(46) 



(47) 



and where ip*i and ip'^ are the values of ip*i and its p- 
derivative at the initial time pi. From Eq. @, ip* is 
expressed in Planck units and we will allow values of ip*i 
to be of order unity. Interestingly, we see that in the 
radiation dominated era, ip* rapidly tends to a constant 
value. The field derivative is just 



V* = 




-(p-Pi) 



(48) 



so that it is divided by e n in Ap = n e-folds. In particular 
if we send pi — ► — 00 then its variation between pi and 
some time in the radiation era is 



Aip* — > — \/31n a 4 + J 1 



(49) 



6 



It follows that, as long as <p' ti <C v3, |A</j*| ~ and the 
field gets frozen at a constant value during radiation era. 
These properties can be recovered easily from the form 
of Eq. 112 1H of the evolution equation since it implies that 
<p* decreases as RZ^- This behavior is quite general for 
dilaton-like fields [43. In conclusion, deep in the radia- 
tion era (much before nucleosynthesis) the initial condi- 
tion can be chosen to be Lp*- m = and ip*i n =constant. 



2. Mass thresholds 

The previous analysis ignores an interesting effect 
that appears when the universe cools below the mass of 
some species \i T ~ m x . This species becomes non- 
relativistic and induces a non-vanishing contribution to 
the r.h.s. of Eq. (|21[1 . For example, during electron- 
positron annihilation, the r.h.s. of Eq. (|21[) depends on 

S e = (p e - 3P e )/p rad . 

In the Jordan frame the total energy density of the 
radiation is 



Slices at constant I 



1 








: / 10 










40 

on — — 






_/ / ou _ _ — — -- r _ — — _ _ 






90 




■If' / / 70 










60/ 






'/ \ 

f; 1 


, 1 ,/, , 1 , , , 1 , , , 1 ,\ 


<!,,,!,, 


1 1 1 1 



0.4 0.6 0.8 



1.2 1.4 1.6 



Prad = 9*{T)—T 4 



(50) 



where g* is the effective number of relativistic degrees of 
freedom, 



9*(T) 



fcrmion 



9u 

bosons 



(51) 



and T is the Jordan frame temperature of the radiation, 
as long as the particles are in thermal equilibrium with 
the radiation bath. The term p e — 3P e takes the form 



-3P fi = 



9e 2 

2^ 



dq 



Intoducing x = E/T and z e 



1 \J q 2 + 
m e /T, we deduce that 



(52) 



£e(T) 
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9*(T)' 



1 



-dx 



(53) 



so that the Klein-Gordon equation (|44|l takes the form 

(54) 



3 - ¥V 3 



The force term depends on the temperature which de- 
pends on A(ip*) and p. When this term is no longer ef- 
fective, the field evolves according to Eq. I|45ll and tends 
to another constant, <p* ou t- The relation between <p*- m 
and v?*out, or equivalently between a ln — a((p^ ln ) and 
aWt j has a complicated structure. Eq. (|44|l is almost 
a damped oscillator (because of the non- linear term). 
When E e (T)/3<^*i n is small, the field does not have the 
time to oscillate while S e is non-negligible, and the rela- 
tion ai n -a ou t is linear. For larger values of (3 and/or a ln 
one gets damped oscillations so that a out < a- m . Figure 
depicts the relation ai n -a ou t for various values of the pa- 
rameter j3 and Figure [21 illustrates the complexity of the 
full solution of this equation. 



FIG. 1: fflout as a function of a[ n for different values of (3 
between 1 and 100. We see that a ou t < aj n which reflects the 
attraction towards general relativity during electron-positron 
annihilation. 



3. Details of the field dynamics near threshold 

Let us now investigate the dynamics of the attrac- 
tion toward GR during a mass threshold in more detail. 
In general, the temperature is related to the integration 
variable p by 



T[ip*,p] = T 



A 



g 7 ( T o) 
97 (T) 



1/3 



(55) 



where g 7 is the effective number of relativistic particles 
entering the definition of the entropy, taking into account 
only particles in equilibrium with the photons. 

The dependence T[ip*,p] and the non-linear term (p* /3 
make Eq. H54|) difficult to integrate. In regimes where j3 
and <£>*; n are not too large then it can safely be approxi- 
mated by 



(56) 



This approximate equation assumes that the field is slow 
rolling and that A(ip*) does not vary much during the 
transition. It is a linear equation in tp* so that its solution 
is proportional to p„ n . 

Fig. El compares the solutions of the two equations (ap- 
proximate and exact) for a single mass threshold. Indeed 
as long as (3 is small, the field is slow rolling and A re- 
mains almost constant during the transition. This is seen 
in the top panel of Fig. for (3 = 1 and the field evolves 
to ip„ « 0.84</?* nl . However when (3 is large and as a con- 
sequence A(ip„i n ) is also large, the variation of ip„ during 
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FIG. 2: The general structure of a out as a function of a lrl and 
j3. This illustrates the complexity of the solutions of Eq. (14411 



the transition implies, because of the relation 1)55(1 that a 
given width, AT, corresponds to a larger Ap, while this 
latter is fixed for the approximate solution. This implies 
that the attraction toward tp* = is more important. 
This progression is seen in the middle and lower panels 
of Fig. El 

Fig. 0| compares the value a ou t/aj n as a function of , 
When Eq. (|56(l is used, we recover the result of Ref. 
In this case, since this equation is linear, a ut/oi n does 
not depend on the initial value of ip*i n and is a universal 
function of (3. This is compared to the ratio obtained 
from the integration of Eq. I|54() . As long as tp*i n is small 
(typically of order 0.1), both results agree (because <p' 
remains small and A does not vary significantly). How- 
ever, aout/fliu is typically 10 times smaller when ip*- ul is 
of order unity. This agrees with the results depicted in 

Fig.m 

In conclusion, we see that both the depedence of the 
source term for ip* and the non-linear term in ip'^ lead to 
significant modifications of the dynamics when the initial 
value of the scalar field or j3 are large. 



4-. Expected value of 

In the previous analysis we have restricted ourselves to 
ai n between and 3, mainly for numerical reasons. We 
can now in a position to justify this choice. Indeed, and 
as noted earlier, because the scalar field is frozen during 
the radiation dominated era, we need only specify <Pn n 
as an initial condition. For a given value of (3, this fixes 
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FIG. 3: Evolution of the scalar field in phase space during a 
transition for /3 = 1, 5, 50 from top to bottom when we assume 
<p*m — 1- The solid line corresponds to the exact solution of 
Eq. 1541 while the dashed line corresponds to the solution of 
the approximate equation I56H . 



the initial value of ai n . 

It is difficult to predict the value of a- ln from general 
arguments. For instance, if we expect 99* *~ 1 (in Planck 
units) at the end of inflation, this means that a- m ~ f3 
and <Zi n ~ {3/2. In this case, one would indeed like to 
investigate a m with values up to roughly 50. On the other 
hand, if we expect a deviation from general relativity of 
order one at the end of inflation, then we might expect 
ai n ~ 1, or (p* r~j P^ 1 and a- m ~ /3 _1 /2. In that latter 
case restricting to a; n ~ 3 would be safe. Clearly, without 
a detailed model of the inflationary period it is difficult 
to determine the "natural" range of variation of aj n . 

To get some insight on the expected order of magni- 
tude of Oi n just before the period of electron-positron an- 
nihilations, we must investigate the effect of higher mass 
thresholds. To that end, we consider an extention of 
Eq. I|54|) in which the source term is replaced by a sum 



S(T)= ]T Ei(T), 



(57) 



species 
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FIG. 4: Evolution of a ut/ciin as a function of /3. (Top): when 
we used the approximate equation 15611 . it does not depend on 
the initial value of the scalar field, ip»i u . (middle and bottom): 
we use the exact equation 1541 . This equation being non- 
linear the ratio depends on the initial value of yj*m. 

where 

^ = ?^ 2 £^ di < 58 » 

with Zi = rrii/T, e = +1 for fermions and e = —1 for 
bosons. 

In principle, all massive standard model particles will 
play a role. In addition to electrons and positrons, 
we must consider the effects of muons, pions, charmed 
quarks, taus, bottom quarks, bosons, Z° boson and 
the top quarks. The role of lighter quarks is tied to the 
quark hadron transition whose effect we do not include. 
Nor do we include the effect of the Higgs boson due to 
its as yet uncertain mass. Fig. [5] depicts the evolution of 
£ with (the Jordan frame) temperature. In particular, 
it shows that the effects of the various thresholds cannot 
be considered separately because the transitions overlap 
and the scalar field, does not have time to settle back 
to tff^. = between two transitions. Also note that, fortu- 
nately, the last threshold (electron-positron annihilation) 



is almost decoupled from the previous ones. Thus, we will 
be able to compute the state of the scalar field just be- 
fore the last transition which is of primary importance 
for BBN. 
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FIG. 5: The source function, E(T), entering the Klein-Gordon 
equation when the mass thresholds corresponding to the par- 
ticles listed in the text. The dashed curves show the individual 
particle contributions, Ej(T), and the solid curve shows the 
sum, E(T). 



Fig. describes the dynamics of this multi-threshold 
phase (including the electron-positron annihilation). As 
long as (3 or tp*i n remain small, we see that each of the 
four peaks of E, corresponds to a well defined departure 
from = with movement towards smaller propor- 
tional to the initial value ip„ n . For larger values, the field 
is first slow-rolling and then oscillates around = 0. In 
this case, we see that the effect of the four peaks cannot 
be considered separately because the field does not have 
time to settle back to tp'^ = between two transitions. 

The evolution of ip*, when mass thresholds are non- 
negligible, allows us to determine the value of aj n prior 
to the period of electron-positron annihilation. We de- 
note this value by a cc . Fig. [7| shows a ee /a- m , that is the 
value of a(tp*) just before electron-positron annihilation 
compared to its initial value at very high temperature, as 
a function of (3. The attraction toward general relativity 
is very drastic. In the case where ip^ m is of order unity, 
we conclude that a cc < 1(T 4 x a in - f CT 4 /3/2 < 5 x f (T 3 . 
It follows that restricting to a; n = 0—3 at before electron- 
annihilation is a safe limit even if ~ 0(1) at the end 
of the inflationary phase. 

Note also that phase transitions are another source of 
attraction toward general relativity. We have not in- 
cluded either the quark-hardon transition or the elec- 
troweak transition in the previous analysis. During a 
phase transition, there is generally a significant modifi- 
cation of the equation of state which will induce a source 
term in the Klein-Gordon equation. 
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FIG. 6: Dynamical evolution of <p* when the source term is 
described by Fig. |S] (Top) we have set /3 = 1 and p^ lrL = 
0.1,0.5,1 (solid, dashed, dotted). (Bottom) We have fixed 
(/3»i n = 1 and taken ft = 1, 15, 20, 50 (thin solid, solid, dashed, 
dotted). 
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FIG. 7: Evolution of a ee /ai n as a function of /3 when the source 
term is described by Fig. [KJjust before the electron-positron 
annihilation. (Top): (3=1 (Bottom): (3 = 0.1. 



5. Radiation-matter transition 

In principle, BBN will place a constraint on the value of 
a OU f As such, our constraint will in effect be dependent 
on ai n which is unknown. To compare these constraints to 
the ones obtained in the Solar system, we need to relate 



a ou t to a,Q. We allow the code to integrate the evolution 
equation up to the present, so that we obtain ao directly. 

For the particular case of a vanishing potential or as 
long as the field is slow rolling, ip' <C 1, one can approx- 
3 so that Eq. (|4*4l takes the simplified 



imate 3 
form 



'2 



dy 2 



(59) 



where we have introduced the variable y = i?*/i?dcc* 
and used the fact that the gas is a mixture of presure- 
less matter and radiation and the equation of state is 
w = 1/3(1 + y). This equation allows us to relate the 
value of the scalar field deep in the radiation era but af- 
ter BBN, ifout*, to its value today, ipo*. Its solution is 
a hypergeometric function, fp(y) = 2-Fi[s, s*, 2; —y] with 
s = 3/4 



iy/3(P-3/8)/2 so that 

fa* = (p<mt*fp(ya) 



(60) 



where the matching to the analytical solution has been 
performed after the end of nucleosynthesis at a time 
where tp* is constant. y is given by 



Vo 



Ro A cq 
-Rdcc Ao 



(l + z eq )exp[(a2 Q -a%)/20\. (61) 



This method avoids integrating the system I|37I41|I to the 

present but requires a determination of yo- Indeed when 
tp has not varied significantly between BBN and equality, 
then 



y ^ (1 + z cq ) exp[( 



«o)/2/3] 



(62) 



However, this solution cannot be generalized to a A-CDM 
or to extended quintessence models. For this reason we 
do not use this method and integrate the system numeri- 
cally from Zi n to z — 0. Figure [S] compares our numerical 
integration, from which we determine the exact value of 
yo and the analytic solution JSUJ). We see that the agree- 
ment is almost perfect. It can be checked that an error 
smaller than 10% on the evaluation of yo left o-o almost 
unchanged. Let us emphasize that in more general cases, 
i.e. for different potentials and coupling functions, such 
an analytic solution is in general not known so that the 
full numerical approach is necessary. 

The solution (|60|) implies that ip'^ = (poutgpiyo) with 
3/3(2/o) = -3/%0 2-Fi(l + s, 1 + s*,3; -y )/4. It is then 
possible to estimate, from Eq. I|32|l . the value of ao as a 
function of (a ou t,/3), 



1 



cr = 2a out g/3(?/o)- 



/3 

(i+«Lt)/J(ao) 



1 



^ff-fp (2/0)5/3(2/0) 



(63) 



As shown by Fig- EI as soon as a ou t < 1, the constraint 
on (T is satisfied. This means that for the quadratic 
coupling model, nearly all parameter choices satisfy this 
constraint. 
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6. Equivalent speed-up factor 




numerical 
analytical 




FIG. 8: Comparison of the numerical integration and the an- 
alytical solution for a fiat CDM model with a ln — 1. 




Log 



As long as V = and one can neglect the curvature 
term, the Friedmann equation l|18|l can be written, using 
Eq. (23> as 



'(i + «^) 2 



H 2 = 8irG*pA 2 



(64) 



Comparing this to the standard Friedmann equation, 
3H^ R = 87rG*Ag(l + al)p, one obtains the speed-up 
factor defined to be the ratio of the Hubble parameters, 



A(<p*) 1 + a(ipM 



1 



.4 



y/1 - p' 2 /3 y/T+C% 



(65) 



Figure shows the variation of the speed-up factor dur- 
ing BBN for various values of (3, taking into account the 
effects of electron-positron annihilation. £ is constant 
above z ~ 2 x 10 11 and below z ~ 10 9 . For large val- 
ues of /3, typically j3 > 5, the attraction toward general 
relativity is so efficient that £ ~ 1 for z < 10°. For 
smaller values, £ is frozen at some constant value £ > 1 
at the end of BBN and will be driven towards 1 only 
when the subsequent evolution due to matter domina- 
tion will be significant. For very small values of f3, as 
pointed out in Ref. 35] and as we have shown earlier, ip* 
is almost constant during the electron-positron annihila- 
tion period. As a result, cp* ut = <P*in and ip'^ ~ so that 



21n£ = /T 1 (a 2 n 



ln(l+a 2 ). 



(66) 



In a more complex situation, one cannot approximate this 
factor by a constant and the full dynamics during BBN 
must be determined. In particular, we see that £ drops 
around the time the neutron-to-proton ratio, n/p, freezes 
out, but generically reaches a constant value during the 
nucleosynthesis period. 

More tuned models in which the variation of £ is not 
finished during BBN may lead to some signatures on the 
primordial abundances (see e.g. [5l| for a proposal). In- 
deed, no model independent statements can be made but 
in general we expect them to be very constrained, in par- 
ticular if the mass thresholds prior to electron-positron 
annihilation are taken into account. Such models can eas- 
ily be discussed in future works with the tool presented 
here. 



FIG. 9: <To as a function of a out and f3. 
contour lines give the value of log \ao\. 



The labels on the 



B. Numerical simulations 



The constraint on ao leading to a ou t < 1 is an a pos- 
teriori argument for not considering very large values of 
aj n at electron-positron annihilation. Very large values 
of a m will in general lead to large values of a out that are 
constrained by local tests on the constancy of the gravi- 
tational constant. 



The time evolution of a((p*) is depicted in Figure ITT1 
for three values of [3. It is obtained by numerically in- 
tegrating the equations of § III Dl by a standard Rungc- 
Kutta method. In the top panel of Figure ^2 The two 
plateaus at high z correspond to the constant values of 
a during the radiation era before and after BBN. Oscil- 
latory behaviour due to the damped oscillation of the 
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FIG. 10: Variation of the speed-up factor during BBN for 
various values of j3 starting from the same value of a- m . For 
small values of j3, £ is almost constant and is driven toward 
1 during the subsequent matter dominated era. For example, 
for p < 0.2 £ is constant during BBN. For larger values of 
f3, the attraction toward general relativity is very efficient 
and £ ~ 1 during BBN. The dashed line represents the time 
of n/p freeze-out and we have indicated the interval during 
which light elements are formed. 



aj\ =0.15, n A = o, p=io 




field, as described by Eq. I|60l) begins when the matter 
density becomes comparable to the radiation density. In 
this case, we find that a ou t «0.05 and ao » 10 ~ 9 . This 
implies that ao ~ 1.4 x 1CP 4 . In addition, d<^*/dp|o = 
ip*/H* w 2 x 10~ 5 leading to a ~ 7 x 10~ 8 , easily sat- 
isfying the post-Newtonian constraints discussed above. 
For the other examples depicted in Fig. ^2 we have 
a out = 1.3 x 1CT 3 , a = 2.7 x IO" 12 , a ~ -2 x IO" 5 , 
ip*/H* ~ 4 x IO" 7 , and cr - -9 x 1CT 10 (/? = 60) 
and a out = 0.5, a = 1.7 x 10~ 7 , a ~ 6 x 10~ 4 , 
ij)*/H* w-lx 10~ 3 , and <j ~ -3 x 10~ 6 (/3 = 1). 

In the middle panel of Figure^J we show the evolution 
of a for j3 — 60. As expected, the larger coupling allows 
for several oscillations during e + e~ annihilation, and sig- 
nificantly more oscillations during the matter dominated 
era. The dashed line shows the redshift corresponding 
to matter-radiation equality. We see that the field starts 
oscillating before equality due to the enhanced coupling. 
Comparing the two panels, we see that as increases 
oscillations begin at higher redshift. 



C. BBN calculations 



£2 m h 2 =0.1S, £2 A = 0, p=l 
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FIG. 11: a(<f>*) as a function of z for ft = 10 (top), (5 = 60 
(middle) and /3 = 1 (bottom) assuming aj n = 1. The dashed 
lines show the redshift corresponding to matter-radiation 
equality. In the top panel, the dotted lines define din, flout 
and ao while in the middle panel it emphasizes the (1 + 2)2 
overall dependence of a(tp). 



The equations displayed in section III Dl have been im- 
plemented in our BBN code [53, ■ The source term 
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due to the electron-positron contribution to the energy 
and entropy density is calculated by a numerical integra- 
tion of the fermi distributions. The calculation starts at 
T=10 12 K well above electron-positron annihilations and 
weak interaction freeze-out with a given value of (3, a- m 
and d<^*/di*=0. For a given value of flbh 2 , a grid of cal- 
culations is performed with a- m ranging from 0. to 3. in 
steps of 10 -3 and j3 ranging from 0.1 to 100. in steps of 
0.1 ((3 < 10) and 1. ((3 > 10). Let us emphasize that 
the range in a- m is conservative given the analysis of the 
previous section. Small steps are needed because of the 
complicated structure displayed in Fig. [21 The 4 He and 
D yields are compared to the allowed intervals, discussed 
in the next section, and for each /3 value, the maximum 
allowed value of a ou t is determined and the numerical 
calculation is extended to obtain the present limits on 
ao and ao- (The minimum allowed a ou t value was found 
to be zero in all considered cases.) However, because of 
the late domination of matter, the limits on ao display 
many more oscillations than the limits on a ou t as a func- 
tion of (3. Hence, steps 100 times smaller in (3 were used 
for the calculation of ao using interpolated values of the 
relatively slowly varying a OU f 

Figures El to El illustrate the dependence of D, 4 He 
and 7 Li in terms of the baryon-to-photon ratio, r\ — 
r]io x 10 -10 , the parameter (3 and the initial condition 
of the field, a m for a model with quadratic coupling and 
vanishing potential. It is evident that the D mass fracion 
is only very weakly affected by the two new parameters 
(f3, ajn). 7 Li is also almost independent of a- ln but the val- 
ley around log(ryi ) = 0.5 becomes deeper as (3 increases. 
The most sensitive of the light elements is 4 He. Its abun- 
dance depends strongly on both parameters. This reflects 
the fact that the expansion rate of the universe is modi- 
fied by this scalar-tensor theory of gravity. 



D. Constraints on primordial abundances 

A comparison of the previous set of computations 
with the observational determination of the light element 
abundances will allow us to set constraints on this class 
of models. The abundance data are obtained from spec- 
troscopic observations and compared directly with BBN 
predictions assuming the WMAP determination of f2b 
lH IHEl IHIll. We now discuss these observations. 



1. D/H 

The best determinations of primordial D/H are 
based on high-resolution spectra in high-redshift, low- 
metallicity quasar absorption systems (QAS), via its 
isotope-shifted Lyman-cv absorption. The five most pre- 
cise observations of deuterium |52j, 

IHHEi in QAS give 
D/H = (2.78 ±0.29) x 10" 5 , where the error is statistical 
only. 




logOV 

FIG. 12: D mass fraction as a function of the baryonic density 
(log(7jio)) and Oi n for (3 = 20. The abundance of D is almost 
insensitive to ai n . 

Using the WMAP value for the baryon densi ty [2q the 
primordial D/H abundance is predicted to be |27j : 

(D/H) p = 2.60±g;i? x 10~ 5 (67) 

This value is in very good agreement with the observa- 
tional one. Nevertheless, as we will see below, the agree- 
ment between predicted D/H abundance and observa- 
tions is not very sensitive to the change the gravitational 
sector of the theory. 

2. 4 He 

4 He is observed in clouds of ionized hydrogen (HII re- 
gions) , the most metal-poor of which are in dwarf galax- 
ies. There is now a large body of data on He and 
CNO in these systems j&g, f° r which an extended 
data set including 89 HII regions obtained Y p — 0.2429 
± 0.0009 H3- However, the recommended value is based 
on the much smaller subset of 7 HII regions, finding Y p 
= 0.2421 ± 0.0021. 

It is important to note that 4 Hc abundance determina- 
tions depend on a number of physical parameters associ- 
ated with the HII region in addition to the overall inten- 
sity of the He emission line. These include, the tempera- 
ture, electron density, optical depth and degree of under- 
lying absorption. A self-consistent analysis may use mul- 
tiple 4 He emission lines to determine the He abundance, 
the electron density and the optical depth. The ques- 
tion of systematic uncertainties was addressed in some 
detail in |5g. It was shown that there exist severe de- 
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FIG. 13: Contour plots for D/H as function of a ln and 7710 
but for different values of /3 : 10, 20, 30 and 40. (For /3= 20, 
this is the same as Fig. I12I 1 The contours are evenly spaced 
by step of Alog(D/H) = 0.2 starting from log(D/H)= -4.8 
on the right. We conclude that the dependence of D mass 
fraction on (3 is mild. 

generacies inherent in the self-consistent method, partic- 
ularly when the effects of underlying absorption are taken 
into account. These degeneracies are markedly apparent 
when the data is analyzed using Monte-Carlo methods 
which generate statistically viable representations of the 
observations. When this is done, not only are the He 
abundances found to be higher, but the uncertainties are 
also found to be significantly larger than in a direct self- 
consistent approach. 

Recently a careful study of the systematic uncertain- 
ties in He, particularly the role of underlying absorption 
has been performed [53 using a subset of the highest 
quality from the data of Izotov and Thuan |56(. All of 
the physical parameters listed above including the He 
abundance were determined self-consistently with Monte 
Carlo methods. The extrapolated 4 Hc abundance was 
determined to be Y p = 0.2495 ± 0.0092 Conserva- 
tively, it would be difficult at this time to exclude any 
value of Y p inside the range 0.232 - 0.258. 

At the WMAP value for 77, the 4 He abundance is pre- 
dicted to be [23 

Y p = 0.2479 ± 0.0004 (68) 

and it is in excellent agreement with the most recent anal- 
ysis of the 4 He abundance [Ejj • As we will show, although 
4 He remains the most discriminatory element for physics 
beyond the standard model, the current large uncertainty 
in its primordial value will impede tight constraints on 
the parameters used to extend minimal Einstein gravity. 



FIG. 14: As in Figure^]but for 4 He. 4 He is the most sensitive 
element to the value of dm- 

3. 7 Li/H 

The systems best suited for Li observations are metal- 
poor halo stars in our Galaxy. Analyses of the abun- 
dances in these stars yields [6(j Li/H| p = (1.23^q'\q) x 
10 |:: . 

The 7 Li abundance based on the WMAP baryon den- 
sity is predicted to be |27j : 

7 Li/H = 4.15±^;|1 x lO" 10 (69) 

This value is in clear contradiction with most estimates 
of the primordial Li abundance, as also shown by |30| who 
find : 

7 Li/H = 4.26l{j;Sj x 10_1 ° ( 7 °) 

In both cases, the 7 Li abundance is a factor of ~ 3 higher 
than the value observed in most halo stars. 

An important source for potential systematic uncer- 
tainty stems from the fact that the Li abundance is not 
directly observed but rather, inferred from an absorp- 
tion line strength and a model stellar atmosphere. Its 
determination depends on a set of physical parameters 
and a model-dependent analysis of a stellar spectrum. 
Among these parameters, are the metallicity character- 
ized by the iron abundance (though this is a small effect), 
the surface gravity which for hot stars can lead to an un- 
derestimate of up to 0.09 dex if \ogg is overestimated 
by 0.5, though this effect is negligible in cooler stars. 
The most important source for error is the surface tem- 
perature. Effective-temperature calibrations for stellar 
atmospheres can differ by up to 150-200 K, with higher 
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FIG. 15: As in Fig. ll3l but for 4 He. Contours correspond from 
left to right to Y p = 0.22, 0.23, 0.24, 0.25 and 0.26 (i.e. black, 
red, green, blue and black respectively.) The abundance of 
He is very sensitive to j3. In addition, the dependence On flin 
is increased for larger values of j3. 




FIG. 16: As in Figure ^] but for 7 Li. As in the case of 
deuterium, 7 Li depends very mildly on a\ n . 



compatible with observations. 



temperatures resulting in estimated Li abundances which 
are higher by ~ 0.08 dex per 100 K. Thus accounting for a 
difference of 0.5 dex between BBN and the observations, 
would require a serious offset of the stellar parameters. 
We note that there has been a recent analysis jH3 which 
does support higher temperatures, and brings the dis- 
crepancy between theory and observations to within 2 
a. 

We are now in a position to directly compare our nu- 
merical results for the BBN production of light elements 
in a scalar-tensor theory of gravity with observations. In 
Fig. El we show the resulting light element abundances 
as a function of a- ln with f2b/i 2 =0.0224 for values of j3 = 
between 5, 10, 15, 20, 25, 30, 50, and 100. Starting with 
D/H, we see from Fig. El that D is always compatible 
with observation as long as (3 > 10. For lower values both 
D and 4 He will set constraints. Of course for very small 
beta, we must have small values of a; n as we approach 
standard GR. In that case, the concordance of D/H is 
also restored. We must also emphasize that 7 Li cannot 
be reconciled with observation in this class of models. 

Fig- El depicts the constraints expected on (a out ,/3). 
The black solid curve shows the maximum possible value 
of a ou t for ai n = — 2. We also show the maximum 
allowed value of a ou t from BBN for two choices of ili,h 2 
= 0.0224 (solid) and 0.024 (dashed) based on 4 He (red) 
and D/H (blue). We see that for all (3, 4 He always sets 
the tightest constraints. Interestingly for j3 > 20, the 
attraction toward general relativity is so efficient that, 
assuming reasonable values for aj n , all abundances are 



E. Using WMAP to fix ft b 

As we have seen from the previous discussion, one can 
set sharp constraints on the primordial abundances if 
flbh 2 is set by the analysis of the CMB anisotropies. It 
is important to note however, that the WMAP data has 
been analyzed in a standard cosmological set up which as- 
sumes general relativity, that is a = (3 = 0. As was shown 
in Ref. the CMB power spectrum in scalar-tensor 
theories is modified in 3 principle ways: (1): the modi- 
fication of the Friedmann equations induces a change in 
the age of the universe and in the sound horizon thus 
shifting the acoustic peak structure, (2) the amplitude of 
Silk damping is modified because it depends on the pho- 
ton diffusion length at recombination and thus on the 
Hubble size at this time, and (3) the thickness of the last 
scattering surface is modified. In any specific model, one 
needs to check to what extent the CMB angular power 
spectrum is modified and decide whether the constraints 
set by WMAP can be used as is or if one needs to go 
through a combined analysis to get new consistent con- 
straints. 

As was shown in Ref. 0] , for the case of the quadratic 
coupling adopted here, CMB anisotropies are not affected 
by this modification to gravity and hence we can safely 
use WMAP data to fix the baryon density. But, in gen- 
eral, this will not be the case in other models (see e.g. 
Ref. [HI). 
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FIG. 17: As in Fig. EDbut f° r 7 Li- The contours are evenly 
spaced by steps of ALi/H=1.10 -10 starting from the value 
of Li/H=1.10- 10 (in black only seen in the lower diagrams). 
7 Li shows little dependence on (3 except for values of log(?7io) 
corresponding to the minimum value of 7 Li. 
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IV. INCLUDING A COSMOLOGICAL 
CONSTANT 



a«L) 



A. Generalities 

The previous model does not account for the observed 
late acceleration of our universe. One can easily gener- 
alize it by introducing a cosmological constant. Let us 
however stress that there is no unique way to introduce 
such a constant in scalar-tensor theories. 

One way to generalize the model is to introduce a cos- 
mological constant in the Einstein frame which corre- 
sponds to a flat potential for the dilaton, so that the 
spin-0 degree of freedom remains massless. In this case, 
we consider models in which 



1 



V((p*) = Vq, a(<p*) = -/?¥>*• 



(71) 



The energy density in the Jordan frame related to the 
constant Vq is not a constant energy density and corre- 
sponds to a potential U(<p) = 2V A- A = 2V F 2 ((p). 

Alternatively, we can introduce a constant energy den- 
sity in the Jordan frame. This amounts to choosing 



1 



V(<p.) = U A 4 {ip*)/2, a{tp*) = -f3ip 



(72) 



FIG. 18: 4 He, D and 7 Li abundances as a function of ai n 
for (3 = 5, 10, 15, 20, 25, 30, 50, and 100 (The maximum 
deviations correspond to the lowest values of 0). The baryonic 
density is set by WMAP observations and is assumed to be 
independent of the scalar component of gravity. For j3 > 
10 D is always compatible with observations and will set no 
constraint while 7 Li cannot be reconciled with observations. 
We conclude that BBN constraints will mainly arise from 4 He. 



The properties of these models can be discussed by 
generali zing Eq. I|44|) when the potential does not van- 
ish [13, llffi As such, we set p v = V/4irG* and P v = 
—p v and p T = p*+Pv, Pt = P*+Pv- Using -0* = 
we obtain 



Pt 
Pt 



-cx(ip* 



3P* 



-ay 



Pv 



Pt 
3P y 



Pt 



(73) 



The value of either Vq or Uq is set by the observed value 
of the cosmological constant density parameter today. 



with ocy = dlnV 1 / 4 /^*. 
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n b h l = 0.0224 (A) and 0.024 (B) 



He constraint 
D constraint 
No BBN constraint 





FIG. 19: Constraints in the plane (a out ,/3) arising from D 
and 4 He BBN abundances. The upper line corresponds to 
the maximum value reached by a out assuming ai n = — 2. 
The middle line corresponds to the constraint obtained when 
D abundances are taken into account while the lower line cor- 
responds the constraint from 4 He abundances. We conclude 
that 4 He alone is sufficient to set the constraints on the model 
and that for /3 > 20, all models with reasonable ai n are com- 
patible with BBN. 



B. Constant potential in the Einstein frame 

In this model, the dilaton remains massless and ay = 
0. It follows that the Klein-Gordon equation for 93* takes 
the same form as Eq. <|44|) . The dynamics of the universe 
is just modified at late time due to the contribution of 
the constant potential to the Friedmann equation (see 
Fig. 1221) . It follows that our previous results for the re- 
lation between a m and a ou t are not affected and that 
the attraction mechanism operates similarly. In particu- 
lar the constraints on the parameters (a ou t,/?) obtained 
in Fig. 1191 remain unchanged and the constraints on the 
parameters (ao, /?) will be modified only by late time dy- 
namics. 

The value of Vq is fixed by 



2Vn 



(74) 



and it dominates only during the last e-fold or so. The 
field is damped during the matter era so that it will be 
slow-rolling close to present time. It follows that tp" <C f' 
so that 

(1 + ClvW* ~ -ao^mat- (75) 
Using WMAP concordance values, one finds <//„ ~ 



FIG. 20: Constraints on (qo,/3) from 4 He observations only 
(See Fig. 1191 using WMAP measurement of Sib- The dashed 
line represents the constraint obtained in the Solar System. 



— 0.2a so that the constraint iFil'l) on the time variation 
of G cav is satisfied, simply because ao is small. 

As a consequence, we expect the effect of an Einstein- 
frame cosmological constant is a shift in the global con- 
straint contour obtained previously when V = 0. Fig- 
ure [21 gives the new bounds set by BBN in the (ao , /?) 
plane. 



C. Constant potential in Jordan frame 

When one adopts a cosmological constant in the Jor- 
dan frame, the dilaton is not massless anymore and the 
value of Uq is fixed by the constraint 



nA0 w 



(76) 



The coupling ay — a. As in the previous case, the po- 
tential will dominate only during the last e-fold. In the 
slow-roll regime 



(1 + fVV*o ~ -C*o(^mat + My)- 



(77) 



With the concordance values, this implies that <p'„ ~ 
— 0.6ao so that, again, the constraint 1M2I) on the time 
variation of G cav is satisfied. As shown in Fig. these 
two models only differ at late times and figure I2T1 gives 
the new bounds set by BBN in the (ao,/?) plane. 
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FIG. 21: Same as Fig. I2UI but with a cosmological constant 
either in Einstein frame or Jordan frame. 
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FIG. 22: Evolution of the scalar field as a function of z in 
models with a vanishing cosmological constant and a cosmo- 
logical constant defined either in Einstein or Jordan frame. 
We see that only the late time dynamics is affected by the 
cosmological constant. 

V. CONCLUSIONS 

This article described the implementation of general 
scalar-tensor theories applicable to a BBN code. The 



formalism allows one to choose any self-interaction po- 
tential as well as any coupling to matter. As such it can 
be applied to models which account for the present day 
acceleration such as (extended) quintessence models. 

The ability to use BBN as a constraint completes 
our set of tools, which include CMB anisotropies and 
SNIa 0| and weak lensing [2]]], to study the cosmologi- 
cal imprints of this set of well-motivated theories of grav- 
ity. All observables are computed using the same formal- 
ism for compatibility. They can be used conjointly to set 
constraints on these theories and on deviations from gen- 
eral relativity during the entire evolution of our universe. 
We emphasize their complementarity since BBN depends 
only on the background evolution and mainly tests the 
attraction mechanism toward GR, CMB is mainly sen- 
sitive to the evolution of the perturbations in the linear 
regime while weak lensing probes the non-linear regime. 

In this article, we have focused on the case of a 
quadratic coupling in order to check our code. In the 
case where V = our results are compatible with pre- 
vious analysis 35] . Note however that they do not rely 
on any specific form of the analytic solution. Also, our 
evaluation of the Fermi integrals that are necessary to es- 
timate the kick during electron-positron annihilation do 
not rely on an approximation but rather on a full nu- 
merical integration. We have used a complete BBN code 
with up to date nuclear reaction rates. Current data 
on the light element abundances have been used to set 
constraints and we have also investigated the effect of a 
cosmological constant on these constraints. For this par- 
ticular model, CMB anisotropies are not affected and we 
are allowed to infer fib from standard CMB analyses. We 
emphasize that in general this has to be checked case by 
case. 

Since our approach is fully numerical, it can be ap- 
plied to any scenario and in particular to extended 
quintessence scenarios such as models with runaway 
fields. In these models, during the radiation era, the field 
evolves to reach a scaling solution. Before this, there may 
be a kinetic phase. According to when this kinetic phase 
ends, various effects on BBN can be expected. In partic- 
ular </j*i n =const. may not be a good approximation. It 
was also proposed that the coupling to dark matter may 
be different to the coupling to standard matter. This 
hypothesis relaxes the Solar system bound and allows 
higher values of a c dm- All these questions, and others, 
will be adressed in following works. 
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